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Abstract 

Supersymmetry and extra dimension need not be mutually exclusive options of physics for the TeV scale 
and beyond. Intrinsically higher dimensional top-down scenarios, e.g. string models, often contain super- 
symmetry at the weak scale. In this paper, we envisage a more phenomenological scenario by embedding 
the 4d constrained minimal supersymmetric standard model in a flat 5d S1/Z2 orbifold, with the inverse 
radius of compactification at the TeV scale. The gauge and Higgs supermultiplets are placed in the bulk. We 
assume that only the third generation matter multiplet accesses the bulk, while the first two generations are 
confined to a brane on an orbifold fixed point. From a 4d perspective, the bulk has N = 2 supersymmetry 
which entails a special non-renormalization theorem giving rise to a significant numerical impact on the 
renormalization group running of various parameters. The brane supersymmetry corresponds to = 1 
which we assume to be broken in an unspecified but phenomenologically acceptable way. Given this set-up, 
we study how the gauge and Yukawa couplings and the = 1 brane supersymmetry breaking soft masses 
run through the energy scale exciting the Kaluza-Klein states at regular interval. In the process, we ensure 
that electroweak symmetry does break radiatively. We confront our low energy parameters with the experi- 
mental measurements or limits of different observables, e.g. LEP lower limits on the lightest Higgs boson 
and the chargino, the {g — 2) of muon, the branching ratio of b ^ 57, and the WMAP probe of relative 
dark matter abundance. We present our results by showing the allowed/disallowed zone in the plane of the 
common scalar mass (mo) and common gaugino mass {M1/2) for both positive and negative 1^1 parameter. 
Our plots are the first 5d versions of the often displayed 4d ?7io-A/i/2 plots, and we provide reasons behind 
the differences between the 4d and 5d plots. 
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1 Introduction 



The yet elusive Higgs boson is the primary target of the CERN Large Hadron Collider (LHC). But, LHC 
is also expected to reveal a new ruler of the tera-electron-volt (TeV) territories. The standard model (SM) 
has so far been remarkably successful in explaining physics up to a few hundred GeV energy scale. But 
theoretical inconsistencies of the SM (like, gauge hierarchy problem) and experimental requirements (like, a 
candidate to account for the dark matter of the universe) suggest that there are good reasons to believe that new 
physics beyond the SM is just around the comer crying out for verification. Among the different possibilities, 
supersymmetry and extra dimension stand out as the two leading candidates for dictating terms in the TeV 
regime. These two apparently distinct classes of scenarios cover a wide variety of more specific models. The 
usual practice from a bottom-up approach is to attach an 'either/or' tag on supersymmetry and extra dimension, 
as if the presence of one excludes the other. A more careful thought would reveal that the relationship between 
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these two is not necessarily mutually exclusive. In fact, the presence of higher dimensions is a common feature 
of any fundamental theory valid at high scale. We will get back to this issue a little later. For the moment, to put 
things into perspective, we recapitulate the chronological evolution of the extra dimensional scenarios without 
invoking supersymmetry a priori. We restrict our discussion to the flat space scenarios, as we are not pursuing 
the warped path in this paper besides mentioning it as an aside. 

Flat extra dimensions were first studied HI in a scenario where gravity propagates in a millimeter (mm) size 
compact space dimension, with the SM particles confined to a 4d brane. The motive was to bring down the 
fundamental Planck scale to about a TeV. Subsequently, it was conceived that the brane where the SM particles 
live may actually have a very small size, like 10~^^ cm ~ TeV~^, leading to the concept of a 'fat brane' 13. 
In the context of the present paper, we stick to the fat brane scenario. What are the experimental bounds on the 
fatness of such a brane, more precisely, on the radius of compactification (R)l For universal extra dimension 
(UED) models 131, in which all the SM particles access the extra dimensional bulk, a safe estimate is R^^ ^ 500 
GeV. More specifically, the g — 2 of the muon [4|, flavor changing neutral cun^ents HH?!, Z ^ bb decay HI, 
the p pai^ameter l3l|9l, and hadron collider studies HOl reveal that ^ 300 GeV. Consideration of 6 — )• 57, 
however, implies a somewhat tighter bound {R~^ ^ 600 GeV 113). Methods to decipher its signals from the 
LHC data have recently been discussed too lT3l . On the other hand, in the non-universal scenario where both 
the SM gauge bosons and the Higgs boson propagate in the bulk but the fermions are confined to a 4d brane 1141 . 
ii^i cannot be below (1-2) TeV (Bl. The reason behind the difference in constraints is the following. The 
KK parity, defined by (—1)" for the nth KK label, is conserved in UED, while it is not a good symmetry in the 
non-universal scenario. As a result, while in the non-universal models KK states can mediate many processes at 
tree level yielding strong constraints, in the UED model, thanks to the KK parity, KK states appear only inside 
a loop leading to milder constraints. In any case, in the presence of supersymmetry, all those analyses need to 
be modified with more parameters, which would expectedly lead to a set of more relaxed bounds on R~^. 

Now, what is the motivation of studying a TeV scale (or, a fat brane) extra dimension scenario! The im- 
plications of such models have been investigated from the perspective of string theory, phenomenology, cos- 
mology/astrophysics and high energy experiments. Such models provide a cosmologically stable dark matter 
candidate lT6l . trigger electroweak symmetry breaking successfully through a composite Higgs ITTl . address 
the fermion mass hierarchy problem from a different point of view ITSlR and stimulate power law renormaliza- 
tion group (RG) running yielding a lower (few tens of a TeV) gauge coupling (neai--)unification scale l20l422fl . 
Besides, the running of neutrino mixing angles generated from effective Majorana mass operator in a 5d set-up 
has been studied both in non-supersymmetric l24l and supersymmetric l25l contexts. 

But, what is the advantage of supersymmetrizing it? We argue that supersymmetry and extra dimension need 
not always be seen as two new physics considered simultaneously. In fact, they may nicely complement each 
other in some situations through mutual requirements. This can be seen as follows: 

1 . From a top-down approach, string theory provides a rationale behind linking supersymmetry and extra 
dimension. The string models are intrinsically extra dimensional, and more often than not contain su- 
persymmetry as an integral part. That said, we must also admit that establishing a rigorous connection 
between a realistic low energy supersymmetric model with string theory is still a long shot, though a lot 
of efforts have already been put in that direction l26l . 

'Generation non-universality in fermion localization imposes > 5000 TeV due to large flavor-changing neutral currents and 
CP violation fT9l . 

^The power law loop corrections are admittedly ultraviolet (UV) cutoff dominated. It has been argued that if the higher dimen- 
sional theory contains a larger gauge symmetry which is perturbatively broken then the difference of gauge couplings of the unbroken 
subgroups is a calculable quantity independent of UV completion 1231 . 
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2. The higher dimensional field theories are non-renormalizable, and computation of quantum corrections 
in such scenarios require the existence of a UV completion. Superstring theories, which contain super- 
symmetry, provide the main hope in this context. In fact, a joint role of 'partial supersymmetry' (broken 
at very high scale) and warped extra dimension/compositeness has been envisaged from AdS/CFT point 
of view to account for both the 'little' and the 'big' hierarchies ||27l . 

3. Even after embedding the SM in an extra dimensional set-up, the scalar potential remains unstable un- 
der quantum correction. Supersymmetrization stabilizes it and ameliorates the hierarchy problem. It is 
interesting to note that by admitting chiral fermions and their scalar partners in the same multiplet tac- 
itly provides a rationale behind treating the Higgs boson as an elementary object. An elementary Higgs 
can be perfectly accommodated in a flat extra dimensional set-up. As a corollary, the upper limit on 
the lightest supersymmetric neutral Higgs is relaxed beyond the 4d upper limit of 135 GeV due to the 
presence of the KK towers of top/stop chiral multiplets, and the hitherto disfavored low tan j3 region can 
be resuiTccted [28]. 

4. Each 4d supersymmetric scenario has its own supersymmetry breaking mechanism. The origin of this 
mechanism may be linked to the existence of extra dimension. In fact, one of the eaiiiest motivations of 
a TeV scale fat brane was to relate the scale of 4d supersymmetry breaking with |[2ll29ll. 

Keeping these in mind, we outline the formalism of a 5d supersymmetric model in an /Z2 orbifold which 
contains the 4d supersymmetric states as zero modes. In section |2j we state our assumptions leading to the 
construction of the 5d model and comment on supersymmetry breaking. Furthermore, we explicitly write 
down the particle content and their 5d Lagrangian and illustrate the KK decompositions of the different 5d 
fields. In section [3l we derive the beta functions of the gauge and Yukawa couplings as well as those of the 
different soft supersymmetry breaking parameters diagram by diagram, pointing out how they are all modified 
from their 4d values due to the presence of KK states. In section |4l we discuss the numerical effects of RG 
running and highlight the reason behind the differences between the 4d and 5d scenarios. We also point out 
under what conditions we can ensure electroweak symmetry breaking. In section [51 we exhibit the numerical 
impact of RG running through plots showing constraints in the ?tio-Mi/2 plane. We standardize our numerical 
codes by reproducing the known 4d plots before encoding the necessary alterations for producing the new plots 
pertaining to 5d supersymmetry. This also enables us to compare and contrast the 4d and 5d allowed/disallowed 
zones. Finally, in section |6l we showcase the essential features we have learnt from this analysis. 

2 5d supersymmetry 
2.1 A snapshot of our model 

We highlight the salient features of supersymmetry in higher dimension and outline below the various assump- 
tions that lead to a calculable phenomenological framework. 

1. We consider a 5d flat space time metric. The 5th dimension is compactified on a S1/Z2 orbifold. Orb- 
ifolding is necessary to reproduce chiral zero mode fermions as a 5d theory is vector-like. 

2. We embed the minimal supersymmetric standai^d model (MSSM) in this higher dimensional set-up (sev- 
eral consequences of such embedding, mainly the effects on gauge and Yukawa couplings' evolution. 
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have been studied in |[20l ). From a 4d point of view this leads to a tower of KK states. The massless 
sector corresponds to the 4d MSSM states. 



3. Since in 5d bulk the fermion representation is vectorial, the two-component spinor Q that generates 4d 
supersymmetry will in 5d be accompanied by its chiral conjugate mirror Q'^. Thus a = 1 supersymme- 
try in 5d con^esponds to two different = 1 supersymmetry, or equivalently, a = 2 supersymmetry 
from a 4d perspective. In fact, all members of a given KK mode must fall into a valid representation of 
N = 2 supersymmetry. In fact, each 4d supermultiplet is augmented by new chiral conjugate states and 
together they form a hypermultiplet. 

4. Here we are talking about a massive representation of supersymmetry, where the supersymmetry pre- 
serving Dirac mass plays the role of central chai^ge for = 2 supersymmetry. This charge is not 
renormalized, as a consequence of which the bulk hypermultiplets do not receive any wave-function 
renormalizatioiH II20II30I . We observe that this N = 2 non-renormalization has serious numerical con- 
sequences in RG evolution of parameters. The most notable effect is the blowing up of the Yukawa 
couplings into the non-perturbative regime around 18 TeV, which we will take to be the cutoff of our the- 
ory. This is below the scale of perturbative gauge coupling unification, which is around 30 TeV. Recall 
that in 5d we encounter power law running which results in early (compared to 4d) unification. 

5. We allow the gauge and the Higgs multiplets access the 5d bulk. Thus far what we said is nothing but a 
supersymmetrization of UED. Only the matter multiplets make the difference. In the UED framework, 
all SM particles access the bulk, and thus even though there are two fixed points, there is no brane. One 
could as well have kept some or all of the fermion generations in a brane at a fixed point; the difference 
would be that the scenario would cease to be universal. In the present supersymmetric context too we 
have the freedom of keeping some or all of the matter multiplets at an orbifold fixed point. We note that 
unless we confine at least two generations of matter multiplets on a brane, the requirement of perturbative 
gauge coupling unification leads to a constraint R^^ > 10^^ GeV [1221 . spoiling its relevance for LHC. 
On the other hand, unless we keep the third family of matter multiplet in the bulk we cannot ensure 
electroweak breaking. In view of the above, we let the third generation matter multiplet access the bulk, 
but fix the first two generations at the y = brane. 

6. A^ = 2 supersymmetry forbids Yukawa interaction in the 5d bulk as this interaction involves odd (three) 
number of chiral multiplets. Therefore, we localize Yukawa interaction at the orbifold fixed point where 
the supersymmetry corresponds to A^ = 1. 

7. Now we come to the important question as how we break the residual A^ = 1 supersymmetry. Different 
ideas have been advanced for its realization. One way is to break it by the Scherk-Schwarz mechanism 
|[3ll in which fermions and bosons satisfy different periodic conditions over the compactified dimension. 
Explicit realizations towards this using a TeV-scale orbifold can be found in |[32l . Another interesting 
approach was to break the residual supersymmetry by a second compactification on an orbifold with 
two reflection symmetries, viz. 1(7,2 x Z'^ |[33l . This can be viewed as a discrete version of the 
Scherk-Schwarz mechanism. Both these scenarios yield soft masses which aie UV insensitive due to 
the non-local nature of supersymmetry breaking. From a completely different viewpoint, supersymmetry 

^In other words, for N = 2 supersymmetry it turns out that mn = ?ns, which is analogous to gpt. = Qb for N — A supersymmetry. 
Here m is the Dirac mass (central charge) and g is gauge coupling, while R and B are labels for renormalized and bare quantities. 
Since the Dirac mass of A*" = 2 hypermultiplets appears on the right-hand side of the anti-commutation relation of the conserved super- 
symmetry charges, this mass cannot be renormalized. This is intertwined with the observation that only those terms are renormalized 
which can be written as integrals over all superspace volume. The kinetic term of A'^ = 2 hypermultiplets cannot be written as any such 
integral (see discussions and related earlier references in 1301 . 
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breaking may be infused from the brane-bulk interface |[34l . or transmitted from a distant brane |[35l . 
or arisen from a gaugino mediation set-up [i36il (possibly with a much lower cutoff than 10^^ GeV), 
or triggered by some completely unknown brane dynamics, for example, by a spurion F-term vacuum 
expectation value (vev). 

In the context of the present analysis, we keep the exact mechanism of the = 1 brane supersymmetry 
breaking unspecified. We assume that the supersymmetry breaking scale is of the order of the inverse 
radius of compactification, for example c/ R, where c is an 0(1) dimensionless parameter. 

8. Our main goal is the following: Just like in the conventional but constrained version of 4d supersymmetry 
one starts with a common scalar and a common gaugino mass at high scale (e.g. the GUT scale) and then 
run them down using the MSSM beta functions to find the weak scale spectrum, we do exactly the same 
here by assuming a common scalar mass (ttiq) and a common gaugino mass (M1/2) low cutoff scale 
(18 TeV) and follow the running using the KKbeta functions through successive KK thresholds to obtain 
the weak scale parameters. By adopting a phenomenological approach, we scan ttiq and over a set 
of values c/R, with c varying in the range 0.1 to 1 and R~^ fixed at 1 TeV. 

2.2 Multiplet Structures 

As mentioned in the Introduction, from a 4d perspective the KK towers of matter and gauge fields rearrange 
in the form of = 2 hypermultiplets. A judicious choice of Z2 parity of the 5d fields allows us to break the 
N = 2 supersymmetry to = 1 supersymmetry. We briefly review below the multiplet structures of the fields 
following the prescription suggested in |[37l . 

2.2.1 Vector hypermultiplet 

The 5d super Yang-Mills theory contains a 5-vector gauge field, a 4-component Dirac gaugino and a real scalar. 
When dimensionally reduced to 4d, the gauge field splits into a 4-vector and a scalar, the gaugino splits into 2 
Majorana gauginos, and we still have the real scalar previously mentioned. All these fit into a vector multiplet 
and a chiral multiplet in = 1 language. If we represent the N = 2 vector supermultiplet by V, the 4-vector 
gauge field by A^, the gauginos by A and "0, and define a complex scalar field = ^(S + iA^), where S is the 
5d real scalar and A5 is the 5th component of the 5-vector field, then one can schematically represent the 5d 
vector supermultiplet as 

From a 4d perspective (where the compactified 5th coordinate y is just a label), and in the A^ = 1 language, 
one can visualize the vector hypermultiplet by a vector multiplet V (first column) and a chiral multiplet in the 
adjoint representation by $ (second column): 

V{x,y) = -9a^'eA^{x,y) + ie^eX{x,y)-i9^9X{x,y) + ^9^9^Dv{x,y), 

= (^{x,y) + V2 9i;{x,y) + 9^F^ix,y). (2) 
The Z2 parity of V is so chosen that the V contains a zero mode, but $ does not have any zero mode. 
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The gauge invariant action may be written as (f (fx = f d^x f dy) 



gauge 



d^X 



d^e + h.c. + / d 



V2 



(3) 



where the VF"(x, y) is the field strength superfield con-esponding to V(x, y). 



2.2.2 Higgs hypermultiplets 

From the = 1 perspective, the N = 2 hypermultiplet splits into two chiral multiplets. Thus have a 
hypermultiplet and a hypermultiplet. We can represent them as (the tilde symbol represents Higgsino) 



(u,d) 



HL(u,d) H^i^u,d) 
HL{u/d) Hfi(^u/d) 



(4) 



If we denote the two chiral multiplets inside the hypermultiplet H{x, y) as h{x, y) in left column and h'^{x, y) 
in right column, then one can expand the chiral superfield as 



h/h' = Hl/r + V2 9Hl,r + 9^F, 



(5) 



We assign even Z2 parity to h so that it has a zero mode, and odd Z2 parity to h'^ which does not have zero 
mode. The free action of the hypermultiplets can be written as 



'Higgs 



d^X 



d^9 ( h h" + hhj + { j d^9 h" {85 + m)h + h.c. 



(6) 



2.2.3 Matter hypermultiplets 

Matters have similar" hypermultiplet structures similar" to Higgs: 

where, Fl = {4>L: iPl) {Z2 even) and Fr = {4>r, iPr) {Z2 odd) represent the two = 1 chiral multiplets. The 
free matter hypermultiplet action will be similar to Eq. Q. There are five matter representations, two SU(2) 
doublets Q and L and three singlets u,d,e, where the symbols have their standard meaning. 



2.2.4 Gauge interactions 

When the hypermultiplets are charged under gauge symmetry, their free action can be promoted to take care of 
the interaction in the following way: 



'1 



int 



d^X 



d'^e [Flc'^Fl + Frc-^Fr) + <! / d'OFL ( m + 55 - ) Fr + h.c 



1 



(8) 



where, V = VT" and <^ = ^"■T'^ are Lie-algebra-valued gauge and matter superfields. 
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2.2.5 Yukawa Interactions 



We denote the Yukawa part of the superpotential by Wy, which contains the usual chiral superfield combina- 
tions QHuU, QH^d and LH^e. Since an Yukawa interaction involves three (i.e. odd number) chiral superfields, 
it is not possible to write a bulk Yukawa interaction maintaining N = 2 supersymmetry. For this reason, we 
confine Yukawa interaction at the branes even though the associated superfields are in the bulk. A generic 
action (e.g. involving the third generation bulk matter) can be written as 

^Yuk = / (/ ^y) [-^(2/) + ^iy - ^R)] ■ (9) 

As the Z2 odd fields vanish at the fixed points, they do not contribute to Yukawa interactions. We note that the 
first two matter generations can be kept either at either of the two fixed points. Also, the above Yukawa action 
could have been confined at either of the two branes instead of symmetrically on both. 



2.3 KK decomposition of fields 

In order to obtain the action in terms of 4d component fields, we need to write down the KK decomposition of 
the 5d fields in terms of zero modes and higher KK modes 131 . Each 5d field is either Z2 even or Z2 odd. Only 
the even fields have zero modes. The decomposition of the Higgs fields will be exactly like the matter fields. 



ny 



n=l 



^ n=l 

T t \ ^ t(o) ( \ ^ 2 („) ny 



(10) 



n=l 

00 



* n=l 



3 RG evolution and derivation of the beta functions 

The technical meaning of RG evolution in a higher dimensional context has been amply clarified in |[20l . and 
we merely reiterate it in the present context. The multiplicity of KK states renders any such higher dimensional 
scenario non-renormalizable. So 'running' of couplings or parameters with the energy scale does not make 
much of a sense. Rather, one can estimate the finite quantum corrections that these couplings/parameters 
receive whose size depends on some explicit cutoff A. The contribution comes from KR number of KK states 
which lie between the scale of the first KK state, which is 1/R, and the cutoff A. With this interpretation 
of RG running, we compute the one loop beta functions of the gauge and Yukawa couplings and various soft 
supersymmetry breaking masses. We make the following observations: 

1 . The contribution to the beta function from a given KK mode does not depend on its KK label. 
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2. When we consider different KK thresholds we neglect their zero mode masses, i.e. we assume that the 
nth level KK state is kicked into life when we cross the energy scale n/R. 

3. As we cross different KK thresholds, the beta functions also change. The beta function of the quantity 
X at an energy scale Q, where n ^ QR < {n + 1), can be written as {t = ln((5/(3o)> where Qo is a 
reference scale, e.g. the electroweak scale) 

dX 

— - = /3x , where (ix = Pox + n/3x ■ (11) 
ot 

Here /3ox is the contribution induced by the zero mode (i.e. ordinary 4d MSSM) states (which may be 
found, for example, in the review |[38l ) and f3x arises from a single KK mode. Eq. (fTTI) is our master 
equation using which we perform a diagram by diagram calculation for the estimation of ^x for various 
couplings and parameters. 



3.1 Gauge couplings and gaugino masses 



The running of the gauge couplings (gi) and gaugino masses (Mj) are controlled by 



9i 



levr^ 



6° + nbi 



167r2 



6? + nbi 



(12) 



For the gauge groups U(l) (which cooresponds to gi = \/b/2>g' , which unifies), SU(2) and SU(3), 
(33/5, 1, -3) lEH, and h = (26/5, 2, -2) (H, respectively. 



3.2 Yukawa and scalar trilinear couplings 

We recall that = 1 non-renormalization relates the beta functions of the Yukawa couplings {Uijk) to the 
anomalous dimension matrices (7*) of the superfields. This theorem implies that logarithmically divergent 
contributions can always be written in terms of wave-function renormalizations. Generically, T/jj^ may be 
written as 

/3,».-=7;2/"^' + 7^?/^"' + 7^?/^^"- (13) 

The Feynman diagrams showing the KK contributions to the wave-function renormalizations of the scalars and 
fermions are displayed in Fig.[T] The contribution from the gauge sector cancels exactly as a consequence of the 
N = 2 non-renormalization theorem mentioned in section |2] Diagrammatically, the origin of this cancellation 
may be traced to a relative sign between the A^- and (/)-propagators - see Eq. Only the brane localized 
Yukawa interactions contribute to the Yukawa evolution. We also keep track of the fact that the Z2 odd fields 
have vanishing wave-functions at the two branes, leaving the even fields alone to contribute to the diagrams 
in Fig. [U Here we have made a tacit assumption that although the Yukawa interaction is brane localized, only 
one KK level (n) states float inside the loop at a time. This is a technical assumption to ensure calculability by 
avoiding KK divergence which would have arisen while summing more than one KK index in a loop calculation. 

To appreciate the numerical impact of the bulk N = 2 non-renormalization, we first write down the conven- 
tional 4d MSSM beta functions (i.e. those coming from zero mode states in the 5d context) which contribute to 
the evolution of the third generation Yukawa couplings ll38l : 

^0 ^ JJL 



„ J, * 2 o 2 2 

^Vt yt + ybyb- - ^92 - te9i 



levrn"^*^^ ' ^"^^ 3^^ ^2 15^ 
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Figure 1: Feynman diagrams showing the KK contributions to the wave-function renormalizations of the zero mode 
and Hu. Similar diagrams for the other fermions and scalars may be drawn analogously. Here A^^ is a generic gauge 
field and cj) is an adjoint scalar, both arising from a vector hypermultiplet, as explained around Eq. (0}. 



/3° 



Vb 
167r2 

167r2 



Qylvb + ytVt + VrVr - yS'I 

9 

^VrVr + Sylvb - 'igl - -gl 



3gi 



7 2 



(14) 



The corresponding KK contributions are given by 

/3/ = /3^(gi^O) (/ = t,6,r), 
where the vanishing gauge contributions is a direct consequence of the bulk N = 2 non-renormalization. 



(15) 



The effects of the above non-renormalization can also be felt in the evolution of the trilinear scalar cou- 
plings. The relevant Feynman diagrams aie displayed in Fig. |2] Again, for illustration, we first write down the 
contributions to the beta functions from the zero mode (i.e. 4d MSSM) states ll38l : 



/3° 



/3° 



1 



167r2 



at [ ISyt yt + Vhyb - yS-s 



3<7i 



13 2 
15 ■'^ 



9i + '^abvlvt 



' 32 26 
+ Vt ( y 51^3 + 6ffiM2 + -^g\Mx 



1 



167r2 



16 7 

Ofe ( ISy^yfe + vlvt + vlVr - —gl - "igl - —g\ ) + 'i-atvlvb + larVlVh 



/32 



14 



+Vh\ ^giMs + 6(7^M2 + —gtMi 



V 3 



15' 



1 



167r2 



9 



Uylvr + Sy^Vb - 3g^ - ^51 + 6af,yfe*yr + Vr + ^gtM_ 



18 



(16) 



As expected, the beta functions of the soft supersymmetry breaking parameters are proportional not only to 
those parameters but to others as well, since any non-renormalization theorem ceases to work when supersym- 
metry is broken. For the computation of /J^^, we need to keep in mind the essence of Eq. (fTSl ). i.e. the absence 
of gauge contributions in while solving the coupled differential equations. However, that part of the gauge 
contributions (proportional to gf) to the trilinear scalar couplings which multiply the gaugino masses (Mj) in 
Eq. ([T6l ) would still remain while computing the KK contribution. All in all. 



/3a, =/3°,(a/5«^0). 



(17) 
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^(0) 



H 



(n) 




/ 



, , ~ (71) I - (n) 

( 



(0) 



f,(0) X 




^ g(o) 



»(o) 



Figure 2: Feynman diagrams showing the KK loop contribution to the evolution of a trilinear scalar coupling. The 
diagrams contributing to other trilinear couplings may be drawn analogously. 



3.3 Scalar masses 



We make three observations regarding the KK contributions to the evolution of scalar masses (see Fig.[3l): 



1 . The two diagrams in the lower row of Fig. [3]depend on the Yukawa coupUngs. Hence, they are important 
only for the third generation matter fields. 

2. Recall that in the evolution of the Yukawa couplings the KK contributions from the gauge field and 
the complex scalar (p exactly cancelled thanks to the bulk N = 2 non-renormalization. However, their 
fermionic superpartners contribute to the scalar mass evolution and those contributions add up instead 
of canceling out. This happens because these contributions yield gaugino masses which are = 1 
supersymmetry breaking parameters and hence the non-renormalization theorem ceases to be applicable. 

3. Each KK state in the two diagrams in the top row of Fig. [3] contributes twice that of the SM because 
of the doubling of the fermions (this factor of 2 is highlighted in bold-face in Eqs. (ITSl) and (fT9l) below, 
and also marked on the two diagrams in top row in Fig. O). However, each KK state at the lower row 
diagrams contributes the same as in the SM because the odd fermion modes vanish at the brane where 
Yukawa interaction is confined. 



The zero mode contributions to the beta functions can be found in II381 . Below we display the KK contri- 
butions to the beta functions only for the third generation scalars in external lines (S = TT[Yim'^_], where the 
index i runs over the third generation matter states which are in the bulk): 
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Figure 3: Feynman diagrams showing the KK contributions to the running of the up-type Higgs mass. Diagrams con- 
tributing to the evolution of the third generation soft scalar masses may be drawn analogously. The multiplication factor 
of 2 has been explained in the text - see point 3 under section [Ol 
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Note that for the first and second generation soft scalar masses only the second diagram in the top row of Fig. [3] 
would contribute. 



4 Special numerical features of RG running in 5d scenario 

In this section, we highlight the special features of RG evolution in the 5d scenario. We also compare and 
contrast them with the 4d features. For all our numerical estimates we have fixed 1/R = 1 TeV. 

4.1 The Gauge and Yukawa couplings 

The power law running of the gauge and Yukawa couplings has been discussed in Il20ll22l for the non- 
supersymmetric scenario and in ll20l for the supersymmetric case. As far as the Higgs multiplets are concerned, 
there is a crucial difference between our model and that considered in EOl . In our scenario there are separate 
up- and down-type Higgs hypermultiplets - see Eq. Inside each hypermultiplet only the left column with 
label (L) is Z2 even and its scalar zero mode receives a vev, whereas the right column with label (R) is pro- 
jected out. In other words, the hypermultiplet Hu contains the vev Vu and, similarly, contains Vd- On the 
other hand, ||20| contains a single hypermultiplet, each column of which has a zero mode, one to be identified 
with the up-type chiral multiplet which contains the vev Vu, and the other to be identified with the down-type 
containing Vd- While our approach constitutes a straightforward generalization of Hu and Hd from chiral mul- 
tiplets to hypermultiplets, the choice made in ll20l requires non-trivial boundary conditions. These two different 
assumptions lead to significant numerical differences. In our approach, the gauge couplings converge to one 
another but actually do not meet at a single point, while in [20] the gauge couplings do meet at a point. The 
difference in the number of KK scalar- excitations makes the difference between the two approaches. 

Indeed, both gauge and Yukawa couplings exhibit power law running due to summation over the KK states as 
one crosses the energy thresholds. As we have mentioned in section |2l keeping the first two matter generations 
confined at the brane ensures that the couplings remain perturbative even with R^^ as low as 1 TeV. Starting 
from their LEP-measured values at the weak scale, as we extrapolate the gauge couplings using the KK beta 
functions we observe that the three couplings approach very close to one another near 32 TeV, but they do not 
actually meet at any point, as mentioned in the previous paragraph. 

A crucial point of immense numerical significance is that on account of the special N = 2 bulk non- 
renormalization, the third generation matter hypermultiplet kept in the bulk does not receive any wave-function 
renormalization from the gauge hypermultiplet, which we have illustrated below Eq. (fT3] ). As an important 
consequence of this the Yukawa coupling blows up to large (non-perturbative) values around A 18 TeV, 
which we therefore take to be the effective cutoff of our theory. 

4.2 The gaugino and scalar masses 

We assume that at the highest scale A = 18 TeV of our theory, i.e. just before the Yukawa couplings blow 
up, all scalar masses unify to mo and all gaugino masses to Mi^2- Our high scale parameters are then 

mo,Mi/2,sgn(/x) and tan /3 = v^/vd. 



12 



The gaugino mass running is governed by the evolution of gauge couplings. Since gauge couplings nearly 
meet around 32 TeV, the gaugino masses tend to converge also at that scale. But in the present context, as 
mentioned before, we forced the gaugino masses to unify at 18 TeV. Recall that in 5d the running is short but 
fast (power law), but in 4d it is long and slow (logarithmic). This leads to a general expectation that starting 
from a given high scale value, the low scale predictions would be similar in 4d and 5d. But since we forcibly 
unified the gaugino masses in our set-up, earlier than otherwise expected, we obtain a somewhat different set 
of low scale values. The gaugino mass scaling in 5d is shown in Fig. |4l while in the inset is displayed the 4d 
running. A rough comparison of the weak scale ratios of the three gaugino masses is the following: 

Ml, M2, M3 ~ (0.4, 0.8, 3.0) X Ml/2 (in 4d) , 

Mi,M2,M3 ~ (0.7,0.8,2.0) X Ml/2 (in 5d). (20) 

If i?-parity remains conserved, the lightest neutralino remains the lightest supersymmetric particle (LSP), only 
that its mass is heavier than what is expected in the standard 4d scenario - see Eq. (l20l) . 

Fig- m shows the running of the soft scalar masses. The large top quark Yukawa coupling continues to play a 
crucial role as in 4d. A rough comparison of the weak scale predictions in 4d and 5d is: 

~ mg + 5.5My2 (in 4d) , 

~ + 3.5M2/2 (in 5d). (21) 

Even for the brane localized scalars, the 5d model predicts slightly higher weak scale masses compared to 4(£|. 
The relatively compressed spectrum in 5d, as expressed through Eqs. (l20l) and (|2TI) . is indicative of a lesser 
fine-tuning compared to the conventional MSSM. 

During power law running we ensure that radiative breaking of electroweak symmetry does happen at the 
desired scal^ Just like in 4d, only rajj^ turns negative while all other scalars remain positive. Again, the 
large top quark Yukawa coupling drives this phase transition. A point to note is that unless we keep the third 
generation matter in the bulk, electroweak symmetry would never break radiatively in our class of models. 
Furthermore, we had to take in 5d a factor 1 to 3 larger (than 4d) value of fi at the cutoff scale. 

5 The mo — M1/2 parameter space 
5.1 Numerical procedure 

For our numerical estimates we go through the following steps: 

1. We scan rriQ and M1/2 through the range [0.1 — 1.0]/R. We fix tan /3 = 10 and take both positive 
and negative values of fi. Our cutoff is 18 TeV, which is the scale where the Yukawa couplings blow 
up to non-perturbative values. We use one loop RG equations as displayed in section |3] To obtain the 
weak scale parameters we follow the logarithmic evolution like in the 4d desert, but now using different 

''when the TV = 1 supersymmetry breaking is induced by brane-localized spurion vev, the structure of the third generation soft 
masses may, in principle, be different from what we assumed due to KK mixing. For calculational simplicity we have ignored KK 
mixing which, we believe, would not affect our main qualitative features. 

' Radiative electroweak symmetry breaking has been discussed in the context of some specific realization of supersymmetry breaking 
in an orbifold (33l[39l[40|. 
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beta functions for different energy intei^vals separated by successive KK thi^esholds - see Eq. ([TT]) . So, 
operationally, power law running is treated by taking many logs. 

2. For each input combination, we perform a consistency check to ensure correct electroweak symmetry 
breaking, and accept only those inputs which admit this phenomenon. 

3. We then feed the weak scale spectrum into the code micrOMEGAS f4Tl, and using this software package 
calculate the dark matter density (JIdm), Br (6 — )• 57), Aa^ = {9 — 2)^/2, and Ap. Since we consider 
1/R = 1 TeV, even the lightest KK states ai^e somewhat heavier than the lighter section of the zero 
mode spectra. As a result, without any significant numerical compromise we neglect the direct loop 
contributions of the virtual KK particleo In other words, the KK effects feed into the calculation of low 
energy spectra via power law running, but after that we rely on the standard 4d computations encoded in 
micrOMEGAS. This approximation is good enough for our purpose. 

4. We compare the predictions of the above observables with their experimental values/constraints, and 
translate the information into the inclusion/exclusion plots (Figs.[6]and|7]l in the tuq — Mi 12 plane. The 
4d plots have been reproduced to serve as a guide to the eyes for capturing the 5d subtleties. We note that 
our 4d plots are in agreement with the ones in the existing literature, e.g. with B3l . 

5.2 Comparison between 4d and 5d models 

We highlight only the major differences between the 4d and 5d models that appear on the mQ-Mij2 plane. 

1 . We assume that i?-parity is conserved. In the 4d scenario the Ughtest neutralino is the most likely candi- 
date for an LSP. In the 5d model the situation is somewhat tricky. Indeed, the 4d LSP is still an LSP here 
which is the zero mode neutralino. Besides, if the. KK parity remains conserved, then the n = 1 mode 
of photon tower, namely 71, and its supersymmetric partner 71 are also stable dark matter candidates. 
However, the KK parity is unlikely to respected by the brane-bulk interaction. In our numerical analysis, 
we have treated the zero mode LSP as the dark matter candidate. 

2. We have taken a Sa range of the five year average of WMAP dark matter density (0.087 < Odm^^ < 
0.138) II44II . We raise a caution here that //KK parity remains conserved and we have two more dark 
matter candidates, as mentioned above, then the edge of the allowed band arising from the lower limit 
of r^DM would be further stretched. Note further that in the 5d case there is a slight broadening of the 
WMAP allowed strip compared to 4d. This happens because of a combined effect of Eqs. (l20l ) and (|2TI ) 
leading to a reduced sensitivity to Mi 12 variation. 

3. The region where the lightest neutralino satisfies the dark matter constraints extends to a higher value 
of Ml/2 in 5d compared to 4d. What matters most in this context is the mass difference between the 
chargino and the lightest neutralino, i.e. (M2 — Mi). It is clear from Eq. (l20l ) that starting from a given 
Afi/2 this difference is smaller in 5d than in 4d. This explains the shift of the band to the right side. 

4. We have taken 2.65 x 10"^ < Br (6 ^ 57) < 4.45 x 10"'' El, and 10.6 x 10~^° < Aa^^^ = 
{g — 2)^/2 < 43.6 x 10~^*^ B6l . There is nothing much to distinguish between 4d and 5d from these 
two observables. 

^Although b s-y, Z — > bb, Ap, Bq-Bq are all finite upto one loop in UED or its supersymmetric version thanks to the KK 
number conservation at the tree level Lagrangian, in scenarios containing brane-bound fermions one encounters a UV sensitivity even 
at one loop 1421 . A sensible UV completion is needed to counter this. See also i23j and footnote 2 of this paper. 
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5. We have not included the direct loop effects of the virtual KK states for any of the weak scale observables. 
For R^^ = 1 TeV or more, for processes like muon anomalous magnetic moment or 6 — 57, such effects 
are numerically negUgible, but only for the Higgs mass it makes a difference. In Figs.[6]and|7J the entire 
region to the left of the line mai^ked with = 114 GeV which looks otherwise disfavored will be 
resurrected if we include the KK loop correction to the Higgs mass |[28l . 

6 Conclusions and Outlook 

We reiterate once again that supersymmetry and extra dimension are interestingly intertwined both from theo- 
retical and phenomenological perspectives. The presence of extra dimensions is an essential part of any high 
scale fundamental theory, and supersymmetry is quite often an integral component of such theories. Supersym- 
metry, on the other hand, provides a sensible UV completion of the extra dimensional models. The connection 
between supersymmetry and extra dimension is further strengthened by the consideration that the origin of 
supersymmetry breaking may reside in extra dimension. In a general class of scenarios in which TeV-size soft 
masses are generated at a relatively low cutoff scale (like 20-30 TeV) with intermediate resonances, one can 
find a common ground where our analysis will be applicable. 

The logarithmic running in 4d from 100 GeV to 10^® GeV is replaced in 5d by fast power law running on a 
shorter interval from 100 GeV to about 30 TeV thanks to the KK states. This is a feature of extra dimension. 
What is technically/operationally special about 5d supersymmetry is that it provides a special N = 2 non- 
renormalization that forces us to consider an early cutoff (~ 18 TeV). 

The constraints in the mo-Mi/2 plane have been placed for the^zr^f time in this paper. The ratio M1/M1/2 is 
higher in 5d compared to 4d. For this reason the allowed region in the 5d plot extends to larger values of Mi /2 
compared to the 4d plot. 

It should be noted that could in principle be much smaller, e.g. a few hundred GeV, than our reference 
value of 1 TeV. In fact, the first KK excitations of SM particles could very well be lighter than zero mode 
superpartners. It is neither our intention, nor it is possible within the scope of the present paper, to distinguish 
the two different options. 

Two issues require further studies beyond the scope of the present paper: (i) Besides the lightest neutralino 
(the usual 4d LSP), ifKK parity remains conserved, there may be two other candidates of dark matter in this 
model. One is 71, the n = I level photon, and the other is its superpartner 71. Although brane-bulk interface 
is likely to frustrate KK parity conservation, yet for the sake of completeness one should study the combined 
effects of all three dark matter candidates, (ii) It will also be interesting to revisit the lower limit on R^^ in a 
supersymmetric scenario, which we suspect would be relaxed. 
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Figure 6: The tjiq — M1/2 parameter space for p > 0. The 4d plot is also shown to guide the eye. We keep tan/3 = 10 
for all the plots. The region ruled out by Br (6 — > sj) is shaded in light green (lightest shade) ,the f LSP region is shaded 
in red (darker shade) and the region favored by (g — 2)^ is the region between the two blue (darkest shade) lines. The 
WMAP allowed region where .087 < fl^Mh^ < -138 is shaded in black. We also show the contours for rrih = WAGeV 
and Tnj^± = 103.3Gel^, the region to the left of these lines are ruled out by LEP exclusion limits. For the 5d models, the 
Higgs contour shown does not include the virtual KK contribution. 
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Figure 7: Same as in Fig.^but for fi <Q. The entire region is now disfavored by [g — 2)^ 
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